RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 5, May 1965 

On Electromagnetic Radiation From a Magnetic 
Dipole With Arbitrary Orientation Embedded 
in a Lossless Magneto-Ionic Medium 

H. Motz 
Department of Engineering Science, University of Oxford, Oxford, England 

(Received October 22, 1964; revised December 7, 1964) 

The methods of an earlier paper for calculating the radiation from an electric dipole are extended to the case of a mag- 
netic dipole, and numerical data concerning the radiation admittance of such a dipole embedded in an infinite magneto-ionic 
medium are given. 

1. Introduction 

In an earlier paper [Kogelnik 1960a] the radiation resistance matrix of an electric dipole with 
arbitrary orientation with respect to the magnetic field was calculated in terms of the tensor ele- 
ments characterizing the medium in which it is embedded. Numerical data concerning the resis- 
tance matrix were published in the Proceedings of the Upsala Conference [Kogelnik 1960b]. The 
extension of the method to the case of a magnetic dipole [Kogelnik, Oxford Thesis, 1960, unpub- 
lished] was given at the Copenhagen Symposium on Antennas [Kogelnik and Motz 1963]. The 
purpose of this paper is to present numerical data which were not given in that paper, and to briefly 
but completely reformulate the theory. 

2. Definitions and Basic Analysis 

We treat a dipole embedded in an infinite, homogeneous nonmagnetic (/x rel = 1) medium. The 
a-c quantities are described by their complex amplitudes. A column vector is expressed by a 
symbol a and its transpose, a row-vector, by a. The Hermitian conjugate of a vector a is denoted 
by a + — a x . Symbols like aa, kk, etc. denote dyades; tensor quantities are denoted by symbols 
a, X, etc. In mks units Maxwell's equations take the form 

VxE = -j4oH-M 

VxH^jWE (1) 

where the introduction of the magnetic current density M will allow us to represent a small wire 
loop immersed in plasma by the equivalent magnetic dipole, with a distribution of magnetic current 

M(r)=jo)p m 8(r) (2) 

where 5(r) is Dirac's 8-function. 

In (1) we have dropped factors e> w expressing the time dependence and e is of the form 
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encountered in magneto-ionic theory. For the source-free (M = 0) medium one easily finds, from 
Maxwell's equation, the wave equation which, when applied to waves of the form E = E exp 
( —7k • r), leads to 

r^(k,k )E=0, (4) 

where 

W E (K ko) = k k - (k k)l + k\l (5) 

For nonvanishing fields one must have 

det W E (K ko) = 0. (6) 

One can easily show that, with e of the form (3) 

det We = k 2 (e, sin 2 6 + e 3 cos 2 6) (k 2 - k 2 ) (k 2 - 2 n ) (7) 

with refractive indices kj/kl, k 2 u \k\ for the ordinary and extraordinary wave given by 

kj u = (gi ~ gj) sin 2 6 + 6163 (1 + cos 2 g) ± [(e 2 - ej - e^f sin 4 6 + 4e 2 e 2 cos 2 g] ^ 2 
k% 2(€i sin 2 ^ + e 3 cos 2 0) 



3. Dyadic Green's Functions and Their Fourier Transforms 

We want to express the radiation field in terms of a given source distribution. Since Maxwell's 
equations are linear we can certainly write 

E(r) = / dr'G E M(r, r')M(r') H(r) - / rfr'&n(r, r')M(r') (9) 

where Gem, Ghm are matrices which may suitably be called dyadic Green's functions. Once they 
are known, the problem is solved. We can obtain them by algebraic operations by means of Fourier 
analysis of the fields of the form 

E(r) = J dkE k e-* kr > (10a) 

E k =8^|drE(r)^ kr) . (10b) 

We also define matrices gEM, gmt, which are Fourier transforms of the dyadic Green's functions 
defined above: 



, = 8^J dk & m 



Mk)e- Mr - r '\ (11) 

Introducing (10a, b) and (11) into (9) we obtain 

E k = | £M M k 

H k =g™M k . (12) 
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4. Computation of the Dyadic Green's Function Transforms 



We shall reduce the computation of the transforms gusi, £em to matrix operations. To do this, 
it is convenient to define a matrix 



K = 



(13) 



0, — k :h k 2 

fa, 0, -ki 

k 2 , k u 
such that vector products may be written as matrix products, e.g., 

kxE = £E. (14) 

With this notation the Fourier-transforms of Maxwell's equations may be written in the form 

-jkE k =-ja)fjLoH k -M k 

-jKH k =jcoe {) eE k 



(15) 



The determinant of the matrix K is clearly zero and therefore K does not possess an inverse. In 
most cases, however e _1 exists. We solve for E k and H k by multiplying equations (15) by l/jiofJLo 
and i/jcoeo respectively and obtain 

H k = ^-yKE k -M k ) E k = (-€- 1 /Jw€ol/«H k . (16) 



Substituting these results into (15) we obtain 

(kk-(k-k-)l + ^€)E k =-jKM k (ke- 1 K + kl)H k =j(oe i) M k , 

since 

ii=kk-(k-k)i. 

Comparing (17) with (12) we find the dyadic Green's function transforms 

A A A 

gEM = —jW] i l K, gHM=jcO€ W„ l 

where We is defined by (5) and Wu is given by 

W H = Ke- l K + kl. 



(17) 



(18) 



(19) 



We shall, in the following, only need gHM> For this reason, and also since the expressions for W^ 1 
have already been given in Kogelnik [1960a], we give below only Wjj 1 . One finds, putting k = A;n 
where n is a unit vector 



^ / _ 1 _ nen/c 4 nn - k 2 k%S + A; 4 J(det e) 



det We 



(20) 
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(20a) 



Su = €ie 3 4- n\(e\ — e|) S21 =jnle 2 e :i + 7iin 2 (ef — ef ) 

512 = — ;Vif€2€3 + rci7i 2 (ef — ef) S22 = ei€ 3 + nf( € i ~~ e D 

513 =jn 2 n 3 e 2 es + rcift 2 eie 3 S23 = ^2^36163 —jnin s e 2 es 

531 =—jn 2 n 3 € 2 €3 4- rci7i 3 €i€ 3 

532 = ^2^36163 +7^1^3^263 

5 33 = 2l^€i€ 8 + (^ + 7lI)(€f-€l). 

5. Power Radiated by a Distribution of Magnetic Current 

The power radiated by a volume element cfo; dy dz of the distribution is given by Poynting's 

dP = - %H x Mdx dy dz. (21) 

Hence the power radiated by the magnetic current system is given by the integral 

P=- h J dr' J drM+(r')G+ M M(r) (22) 

over the magnetic current distribution which, in terms of the dyadic Green's function, transforms 
can be shown to reduce to 

P = - 4tt* [ dk(M+g HM M k ) x (23) 

by the application of well-known 8-function identities. 

6. Radiation Admittance of a Current Loop 

As mentioned in the first section, the magnetic dipole is represented by a magnetic current 
density distribution 

M(r) =jo>p m 8(r). 

This corresponds to an oscillating current in a small loop such that the dipole moment p m is given 
by p m = ijlJcPtt, I = Iq sin o)t where a is the radius of the loop. This has an effect different from 
that produced by a rotating electron in its orbit, which is equivalent to two linearly polarized oscil- 
lating currents with a phase difference of tt/2. We have shown in Kogelnik [1960b] that the power 
radiated by two such linearly polarized currents is, in a gyroelectric medium, not equivalent to 
the sum of the power of the separately excited individual antennas (see appendix). 
The spatial Fourier-transform of the magnetic current density is given by 

M, =i^T (24) 

By application of (23) we find that the complex conjugate of the power can be written as a bilinear 
form 

P* = yPm£nPm (25) 
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with 



Zm 



= -^jdkg IIM . 



(26) 



The elements of z m are measured in ohm ! m 2 (admittance-units per unit area.) The matrix z n 
is split into its Hermitian and its anti-Hermitian part 



and the real power radiated is 



Zm r m ~r ]%n 



p — _ + 

* r ^ Pm r mPn 



(27) 



(28) 



We can substitute in (26) for gnu from (18) with the result 

. -joj€ tt (*» . Q , a f« _ f 2 * ,, A- 2 [nen^nA-yc 2 yc^ + ^J(det e>] 

So far the theory is perfectly general and would also apply to lossy media. We now specialize 
the theory to apply to a lossless medium by assuming that 6i, €2, and e 3 are real. Our numerical 
results apply to this case only. The matrices M m , N are denned by 



N 



1 f 2n 
=— nhd(f)'- 

7T Jo 



sin 2 0; 


0; 





0; 


sin 2 0; 





0; 


0; 


2 cos 2 e 



(30) 



M m = M+ = 



Sd<f> 



2e 1 e :i + (e 2 - e|) sin 2 0; - 2/e 2 € 3 cos 2 ; 

3/6*3 cos 2 ; 2e,e :! + (e'f - e 2 ) sin 2 0; 

0; 







0; 4€,€ :f cos 2 + 2(6? - e 2 ) sin 2 



(31) 



Having performed the integrations with respect to <f> and A:, the final results can be written in the 
form 

_ duF 
2Z ( 



r m '- 



j^f 1 



+ (e 3 -e l )u 2 ](kz-k* l ) 



(32) 
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with Zo = {/jLo/eoY 12 , X() — "T - and F m being the Hermitian part 
k$ 

F m = Her m {n€n(jfc? - k$N - k%\k\ - k 3 u )M m + 2k±{k x - A: n )l(det e) } . 

Thus the matrix f m has the form 

r >n], J r m2, 

0, 0, r m3 



(33) 



(34) 
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with the elements 



Fmi 



— —J 1 ' 

2ZqK Jo i 



dufmiju) 



k^e 1 ^(e 3 -€ 1 )u 2 W l -k] l ) (35) 

where ^ = cos 
/mi = Re {fien^f-^) sin 2 0- *ftAf-*g) [26x63 + (e?-ef) sin 2 0] + 2*J(*, - *„) det 1} 

/m2 = Re {-2^(^-^)6263 COS 2 } 

/ m3 = Re {2nln(A:f-A: I 5 I ) cos 2 — 2Ag(Af — Af, ) [2eie 3 cos 2 + (6 2 -6 2 ) sin 2 0] + afcjfo - k u ) det e } . 

7. Numerical Data Concerning the Elements of the Admittance Matrix 

Since the expressions (35) for the admittance matrix are rather cumbersome, we have evalu- 
ated them by means of a digital computer (Ferranti Mercury). Table 1 gives the element r m \. 
Table 2 records values for r m2 , while r m3 is given in table 3. The entries in these tables must be 

Table 1 





62 


• 
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1 \. 


0.1 


0.2 


0.5 


0.8 


0.95 


1.00 




0.001 


0.247 


0.236 


0.165 


0.057 


0.011 


0.003 




.003 


.249 


.238 


.167 


.061 


.015 


.006 




.01 


.254 


.244 


.175 


.072 


.027 


.019 




.03 


.269 


.260 


.197 


.102 


.058 


.051 




.1 


.322 


.315 


.264 


.188 


.155 


.149 




.3 


.473 


.469 


.438 


.395 


.382 


.383 




.5 


.624 


.621 


.603 


.583 


.583 


.588 


1.0 


1.000 


1.000 


1.003 


1.019 


1.043 


1.057 




1.8 


1.601 


1.604 


1.627 


1.680 


1.727 


1.750 




3 


2.502 


2.508 


2.549 


2.637 


2.708 


2.739 




10 


7.754 


7.767 


7.856 


8.032 


8.162 


8.216 




30 


22.756 


22.775 


22.909 


23.168 


23.355 


23.429 




100 


75.259 


75.284 


75.467 


75.817 


76.064 


76.161 




300 


225.26 


225.29 


225.52 


225.95 


226.25 


226.379 




1000 


750.26 


750.30 


750.58 


751.10 


751.46 


751.608 




3000 


2250.3 


2250.3 


2250.6 


2251.2 


2251.6 


2251.8 




10000 


7499.4 


7499.4 


7499.8 


7500.5 


7501.0 


7501.1 



Table 2 





€2 














e 3 


V €1 
-1 ^s. 


0.1 


0.2 


0.5 


0.8 


0.95 


1.00 




0.001 


0.985 


0.941 


0.650 


0.217 


0.031 


0.000 




.003 


.985 


.941 


.650 


.217 


.031 


.000 




.01 


.985 


.941 


.651 


.219 


.033 


.001 




.03 


.985 


.941 


.654 


.225 


.038 


.005 




.1 


.986 


.943 


.663 


.243 


.058 


.023 




.3 


.987 


.947 


.686 


.289 


.113 


.078 




.5 


.987 


.950 


.702 


.329 


.161 


.129 


1.0 


0.989 


0.955 


0.737 


0.408 


0.264 


0.238 




1.8 


0.990 


0.962 


0.777 


0.505 


0.393 


0.377 




3 


.992 


.970 


.822 


.613 


.539 


.537 




10 


.998 


.992 


.962 


.960 


1.017 


1.060 




30 


1.005 


1.018 


1.124 


1.370 


1.588 


1.691 




100 


1.013 


1.051 


1.327 


1.884 


2.310 


2.490 




300 


1.021 


1.082 


1.523 


2.386 


3.017 


3.272 




1000 


1.029 


1.118 


1.745 


2.953 


3.815 


4.156 




3000 


1.038 


1.151 


1.949 


3.476 


4.553 


4.974 




10000 


1.047 


1.187 


2.174 


4.052 


5.365 


5.874 




100000 


1.064 


1.258 


2.621 


5.195 


6.977 


7.660 
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Table 3 





Ca 
















1 ^\ 


0.1 


0.2 


.05 


0.8 


0.95 


1.00 




0.001 


0.001 


0.001 


0.003 


0.004 


0.004 


0.003 




.003 


.001 


.003 


.006 


.008 


.007 


.006 




.01 


.003 


.007 


.016 


.021 


.020 


.018 




.03 


.008 


.016 


.038 


.051 


.051 


.048 




.1 


.019 


.038 


.090 


.127 


.133 


.131 




.3 


.039 


.077 


.184 


.269 


.293 


.295 




.5 


.052 


.103 


.249 


.369 


.408 


.415 


1 


0.075 


0.149 


0.363 


0.546 


0.615 


0.630 




1.8 


0.099 


0.198 


0.483 


0.735 


0.837 


0.864 




3 


.123 


.246 


.604 


.929 


1.066 


1.104 




10 


.192 


.383 


.946 


1.474 


1.714 


1.785 




30 


.264 


.527 


1.306 


2.050 


2.401 


2.509 




100 


.349 


.697 


1.732 


2.734 


3.212 


3.363 




300 


.429 


.858 


2.134 


3.378 


3.977 


4.169 




1000 


.519 


1.036 


2.581 


4.094 


4.827 


5.064 




3000 


.601 


1.201 


2.991 


4.751 


5.607 


5.885 




10000 


.691 


1.381 


3.442 


5.472 


6.464 


6.787 




1000(K) 


.869 


1.738 


4.335 


6.902 


8.161 


8.574 



multiplied by 



2tt 



3Z ( )A_5 



e 3/2 mhos/m 2 



where Z () CL 120 77 ohm, and Ao is the free space wavelength. The values so obtained are expressed 
in mhos/m~ 2 . We also give graphs of these quantities (figs. 1, 2, 3, 4, 5, 6). 1 
The power can be computed from the formula 

■ r 9 Pw r >" P'"- 





Figure 1. Matrix element r mi as a function of— — 1 Figure 2. Matrix element r mi as a function of — — 1 



with 62/61 as a parameter. 



in the range where it does not depend on the ratio 
62/61. 
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Figure 3. Matrix element r m 2 of the admittance matrix 

as a function of— — 1 with 62/61 as a parameter. 
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Figure 5. Matrix element r m'5 of the admittance as a 

63 
function of 1 with 62/61 as a parameter. 
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Figure 4. Matrix element r m2 as a function of ——I 
with 62/61 as a parameter. 



Figure 6. Matrix element of the admittance as a func- 
tion of— — 1 with 62/61 as a parameter. 



The z-axis (3-axis) was chosen throughout as the direction of the applied magnetic field so that in 
the special case of a dipole oriented parallel to the magnetic field, p is given by 

P = /0\p + = (0, 0, Pz ) 
and the expression for the power reduces to 



8. Conclusion 

The data of this paper should be used with caution. When the refractive index of the cold 
plasma theory tends to infinity, the wavelength tends to zero, thus invalidating the assumption 
that the current loop is small. The reactive power, in any case, diverges. It has no meaning 
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r x 


-jr y 





Fy 


r x 











r z 



for an elementary dipole or for any 8-function distribution of current over a conducting surface. 

When there is no propagation, the radiation admittance is zero. 

The author thanks Mrs. R. Baumberg for carrying out the computing and Mr. E. Gilbert for 
checking results. 

The research reported in this document has been sponsored in part by the United States Air 
Force, under Grant 63-33 and monitored through the European Office, ARDC. 

9. Appendix 1 

It was shown in Kogelnik [1960a, 1960b] that the real power radiated by an electric dipole with 
moment p is given by an expression 



i-f* 



Thus for the special case of dipoles which are linearly polarized in directions parallel to x or to y 
the power is given by 

r. r = — r x \Pxr or ^y = ~2~ r v\Pv\ respectively. 

On the other hand, for the case of a circularly polarized moment with 

Pz = y ±jp x =Py 
we obtain 

Pc = <*> 2 \px\ 2 (rx±r y ) 

and this is not the sum of the powers radiated by two aerials when each is excited in the absence of 
the other. The radiation of the aerials shows 'interference' caused by the coupling due to the 
medium. 

10. References 

Kogelnik, H. (1960a), On electromagnetic radiation in magneto-ionic media, J. Res. NBS 64D (Radio Prop.), No. 5, 515-523. 
Kogelnik, H. (Aug. 1960b), The radiation resistance of an elementary dipole in anisotropic plasmas, Fourth Intern. Conf. on 

Ionization Phenomena in Gases, Uppsala, Sweden (North-Holland Publishing Co., Amsterdam). 
Kogelnik, H., and H. Motz (1963), Electromagnetic radiation from sources embedded in an infinite anisotropic medium and 

the significance of the Poynting vector, Proc. Symposium on Electromagnetic Theory and Antenna, Copenhagen, 1962, 

pp. 477-493, ed. E. C. Jordan (Pergamon Press, Oxford, Eng.). 



(Paper 69D5-502) 



679 



